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Abstract 

We study the productions of first radial excited states Dq{2S) {q = u, d, s) and Bq{2S) in 
exclusive semi-leptonic Bqi (q' = u, d, s, c) decays by the improved Bethe-Salpeter method. These 
25" states can be detected through their strong decays to ground mesons, where the strong decays are 
calculated by the low energy approximation as well as the impulse approximation. Some channels 
have ratios of order 10"^: Br{B+ D°{2S)i+ue) x Br{D^{2S) D*tt) ^ (4.9 ± 4.0) x 10"'', 
Br{B^ D-{2S)l+Ui,) x Br{D-{2S) D*tt) « (4.4 ± 3.4) x lO^"^. These channels could be 
measured by the current B-factories. For Ds{2S), we also obtain a relative large ratio: Br{B^ — >■ 
D-{2S)i+ut) X Br{D-{2S) D*K) ^ (9.9 ± 14.9) x 10"^. Although there are discrepancies 
of the full decay width between the theoretical predictions of D^{2S) and experimental results of 
i:>(2550)°, the new detected state £'(2550)° is very likely the D^{2S) state. 

Keywords: Radial excited Dq{2S) and Bq{2S) states; Semi-leptonic decay; Strong decay; Bethe- 
Salpeter method. 
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I. INTRODUCTION 



In recent years, great progress has been made in hadronic mass spectra. There are 
many new states that have been observed, e.g., the new particles D*o(2317) Dsi(2700 
D,j(2860) Q, Vc B, ^(3872) Q, X(3940) J, r(3940) Q, Z(3930) Q and F(4260) ^ 
Some of these new states are P-wave {L = 1) state candidates, such as D*q{2317), and some 
of them are first radial excited 2S {L = 0, n = 2) state candidates, e.g., Dsi(2700) and r]'^. 
Recently, other new 25* state candidates 1^(2550)° and Z}*(2600)° are observed in inclusive 
e~^e~ collisions near ^/s = 10.58 GeV j^. 

Besides the progress in experiment, there are many approaches to study the heavy ex- 
cited states in theory, e.g., the authors of [lo| considered Be decays to excited 2P and 35* 
charmonium states with the relativistic quark model; the authors of calculated the de- 
cay of Be — Xccli^i where X^s was an excited charmonium state with the light-cone QCD 
sum rules approach; the authors of [l^ used generalized factorization together with SU{3)f 
symmetry to predict the branching ratios of Bg — ?■ Mcc + L where Mcs was a charmonium 
state and L was a light scalar; using the ISGW2 quark model, and the authors of [13] studied 
Be — XccM decays, where X^s was a radial excited charmonium t]c{2S) or '4j(2S). 

Although several papers considered the topics of radial excited 25" states 14j-ll7l|. there 
is still lack of theoretical investigation for the radial excited states Dq{2S) or Bq{2S), where 
g is a light quark. One may also note that there is no other heavy-light 25* state, which has 
ever been confirmed by the experiment except charmonium and bottomonium, this means 
we have little knowledge about heavy-light 25* states, so the study of heavy-light 25* states 
will enlarge our knowledge of bound states and QCD. 

There are many methods to detect the heavy-light 25" states experimentally. For example, 
by analyzing the inclusive productions of D^n^, D^ti^ and D*^7i^ systems, the Babar 
collaboration found new mesons like D(2550)° and D*(2600)° Since there are a large 
number of B data in two B-factories, and the LHC will produce large data of Bg and B^, 
there will be a best and convenient way to detect heavy-light 25 states by B^r exclusive 
decays. In theory, the properties of the mesons were studied by many approaches together 
with the Dyson-Schwinger(DS) equation of QCD or the Bethe-Salpeter(BS) equation or both 



of them 



18 



20[ |. In this paper, we will study the productions of heavy-light 25" states in 



exclusive semi-leptonic decays of i?„/ mesons by the instantaneous approximate BS method 
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22|. 



Knowing mass and life time (or full width) is helpful to detect the resonance experimen- 
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TABLE I: The total decay widths (in unit of MeV) of i:>(2550)° treated as the D{2Sf state; 
our results are estimated by the low energy approximation and the impulse approximation (in 
parentheses). 



Ex [9] this Paper [23] [24] [25] [26] [27] 

130 ± 12 ± 13 43 ±23(47 ± 17) 8.4 8 22.1 124.1 45.35 



tally. In a previous letter |23|, we have studied mass and strong decay of some Dq{2S) and 
Bq(2S) states by the BS method. Unfortunately, there are large mass and width discrepan- 
cies between our predicted result of D^{2S) and that of Z)(2550)°, which is the candidate of 
D^(2S) in the experiment jol. Our predicted mass is 2.39 GeV of D^(2S) [23], while Babar's 
result is 2539.4 ± 4.5 ± 6.8 MeV of £'(2550)° Q. The total strong decay width of our calcu- 
lation is 8.4 MeV (note that we missed a parameter 0.5 in all the strong decays whose final 
state involve 7r°), which is much smaller than the experimental value 130 ± 12 ± 13 MeV. 



There are several theoretical approaches to study the strong decays of this new state [23|-|28]. 
We show the theoretical estimated full decay widths and the experimental data in TABLED 
One can see that there are large discrepancies between experimental and theoretical results 



except the result of [26] 



From our previous calculations, we find that there are three main reasons that result in 
discrepancies. First, we have chosen a simple potential-the Cornell potential-in order to 
reduce the difficulty of solving the BS equation, which is very complicated in this work. 
Second, we chose a group of old input parameters. Since there was no information of the 25* 



state mass in the previous letter 23|, we obtained the masses of 25" states by fitting data of 
ground states with old parameters: mi, = 5.224 GeV, nic = 1.7553 GeV, rru = 0.487 GeV, 



TUu — nid = 0.3 GeV, and other parameters that character the potential [29|. Recently, by 
fitting data of charmonia and bottomonia which include higher excited states, and combining 
with the results of decay constants, we give new set of parameters: rrib = 4.96 GeV, rric = 1.62 
GeV, nis = 0.5 GeV, niu — rrid = 0.3 GeV [21]. Third, our results are model dependent and 
we only consider two OZI-allowed channels to estimate the full decay width. Furthermore, 
the theoretical prediction of decay width is very sensitive to the mass of 25* state. 

In this paper, we focus on the decay and production of the 25" states, not on the mass 
spectra, so we can vary the free parameter Vq to obtain the new mass spectra in TABLE |TT] 
and the numerical values of wavefunctions, which are used to calculate the transition matrix 



3 



elements in appendix B. By varying all the input parameters simultaneously within 5% of 
central values, we obtain the uncertainties of masses in TABLE HTl 

Although we focus on the production of heavy-light 2S states in the semi-leptonic decays 
of Bqi, we would like to re-calculate their strong decays by the re-predicted mass spectra. 
The strong decay widths are very sensitive to the kinematic range, so some new strong decay 
channels with higher mass of 25* states can exist, e.g., there is a strong decay with a P-wave 
state involved in final states. Finally, we calculate the ratios of strong decays to reduce the 
effect of model dependence, and estimate the production ratios of Bqi semi-leptonic decay 
to the first radial excited states, which are reconstructed by the ground particles in terms 
of strong decay. 

The rest of the paper is organized as follows. In section II, we show the formulations 
of semi-leptonic and strong decays. We give the hadronic matrix elements of semi-leptonic 
and strong decays in section III. The results and discussions are given in section IV. In the 
appendices we introduce BS equation and give some necessary formulas for the calculations 
of semi-leptonic and strong decays. 

II. THE FORMULATIONS OF SEMI-LEPTONIC DECAYS AND STRONG DE- 
CAYS 

In this section, we present the formulations of Bqi mesons semi-leptonic decay to 25* 
mesons and the strong decays of 25* mesons. 

A. Semi-leptonic decay of Bqi 

As we mentioned previously, the best way to study 25* state is through the exclusive 
semi-leptonic decay of initial heavy meson (5°, 5+, 5° or Be). Here, we take the B^ — >■ 
D~ {2S)i~^i>i (see figured]) as an example to show the formulation. The amplitude of this 
process is 



where Vet is the CKM matrix element, = — is the charged weak current, and P 
and Pf are the momentum of the initial meson B^ and the final state D~{2S), respectively. 
The hadronic part can be written as 




(1) 



{P))=U{P + Pf), + f.{P-Pf) 
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D-{2S),Pf 



FIG. 1: Feynman diagram of the semi-leptonic decay i?^ — ?■ D (2S)l^ui. 



{D^siPfM,\B'{P)) = 0, 



(2) 



where /+, /_ are the Lorentz invariant form factors. 

We define x = E(,/M, y = {P - Pff/M^, where Ei is the energy of the final charge 
lepton and M is the mass of initial meson. The differential width of the decay can be reduced 
to 



dxdy 



feci 



647r3 



2 \ 



/ Mf ^ , \ mj Mi , 

4(2a:(l-^ + ,)-4x^-,j+^(^8x + 4^-3,-^J 



where Mf, are the masses of the meson and the lepton in final, respectively. 

=/?_+ = /+/_,/?__ = /!. 



(3) 



B. Strong decay of 2S mesons 



As an example, we consider the OZI-allowed strong decay D^(2S) — t- D*^7i^ (see figure 
|2]). In this work, we take the instantaneous approximation for the interaction kernel in 
meson; it is fit to describe the double heavy mesons {Bc,ric) and heavy-light mesons {Dg, 
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3l|, but it is inapplicable to double light meson {K^ vr), which have complicated 



internal structure. In this work, we take the reduction formula, PCAC relation and low 



energy theorem to deal with the strong decay as we did in 32| . The strong decay amplitude 
of figure |2] can be written as 

T 



23 



32| 



{D*%Pf^)\u^,^,d\D,siPf)), 



(4) 
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FIG. 2: Feynman diagram of strong decay D {2S) — ?• D*^tt (low-energy approximation) 

d 



D2s{Pf)\ 




FIG. 3: Feynman diagram of strong decay D {2S) — ?■ D*"ir (impulse approximation). 

where Pf, Pf^ and P/j are the momenta of the D~{2S), final states Z)*° and vr", respectively, 
fpj:^ is the decay constant of tt^; we call this method the low energy approximation in this 
paper. 

As a comparison, we also calculate the strong decays using an alternative method called 
the impulse approximation jssl . 34 1. According to this method, the decay amplitude of 
D~{2S) — )■ Z)*°7r~ can be written as (see Fig. [3]) 



T 



{2ny 



Tr {sM)VpAqf)S2{p'2)r]pMS[{p'i)TAPf^^^^^ - Qf - {a'^Pf, - Qfi))} , (5) 

where g/, qji are the relative momentum of quark-anti-quark in D~{2S) and Z)*°, 5*1 (p'^), 
S2{p2) and S[{pi) are propagators, rip^{qf), fjPf-^{qfi) are the heavy-meson BS amplitudes, 
respectively; p'^ = a[Pf + qf, p'^ = a'^Pj - qf, p'[ = a'/P/i + g/i; a[ -- 



nig 



a' 



mu+rric 



and a'n 



rriu+mc ' 



r^(Pj2) is the BS amplitude of n. 
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After instantaneous approximation, equation ([5]) can be written as 

where qf± = (0,g/), qfi± = qf± + ^^+"^^ -P/i±, ^p^{qf±) is the positive energy wavefunc- 

tion for D^s^Pf)- ip^^{q..^j_) is the positive energy wavefunction for D*^. TT,{qf±; Pf2) = 

fi I I I I 



^75 j^i?7r(Q'/x); one can find detailed calculation of BT,{q'jj_) in |33N35|. 

There are two channels for the D-{2S) meson: D-{2S) D*^^- (0^(25) l^O") and 
D-{2S) Z)*(2400)°7r- (0-(25) O+Q-). The strong decay amplitudes can be described 
by the strong couphng constants 

T(D-(25) -> D*(2400)°7r-) = ^^-(25)550., (7) 

where Gj^-i^2S)D*'^-K and Gj^-(^2S)Df-K a-^e the strong coupling constants, and e is the polariza- 
tion vector of meson D*^ . 

With equation ([7]), we obtain the decay widths 

rD-(25)D*o. = P^,Y.\nD-{2S) ^ 5*°"-^'^ 



8nMf^ 

rz,-(25)D5(24oo)o.- = ^|T(Z}-(25) -> 5S(2400)°7r-)p. (8) 

From equations ([1]) and (jl]), we find that the main task of semi-leptonic and strong decay 
is to calculate the amplitudes {D2siPf)\Ji,\B°{P)) and (L'*°(P/i)|u7^75rf|L'2s(^/))- 



III. THE HADRONIC MATRIX ELEMENTS OF SEMI-LEPTONIC DECAY AND 
STRONG DECAY 

A. Hadronic matrix element of semi-leptonic decay 

The calculation of the hadronic matrix element are different from model to model. In 
this paper, we combine the BS method that is based on the relativistic BS equation with the 



Mandelstam formalism 36| and relativistic wavefunctions to calculate the hadronic matrix 
element. The numerical values of wavefunctions have been obtained by solving the full 
Salpeter equation that we have introduced in the appendices. As an example, we consider 
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the semi-leptonic decay ^ D 
matrix element can be written as 

{D:^s{Pf)\J,\B\P)) = 



2S)i~^i'e. In this way, at the leading order, the hadronic 



22| 



dq 



-Tr 



^;;(c.)^^r(^")7.(i-75) 



(9) 



where q (q^) is the relative three- momentum between the quark and the anti-quark in the 
initial (final) meson and qf = q — ci\T. M is the mass of 5°, r is the three-dimensional mo- 
mentum of D'{2S), {pp~^{q) is the positive Salpeter wavefunction of 5° meson and ^Pp^{qf) 
is the positive Salpeter wavefunction of D^{2S) meson, ip^~^ = 7o(<y5^^'*")^7o. We have given 
the Salpeter wavefunctions for the different mesons and form factors in appendix B. 

B. Hadronic matrix element of strong decays 

We have obtained the amplitude of strong decay D^(2S) Z)*°7r~ in equation (jlj), and 
the key factor is to calculate the hadronic matrix element {D*^{Pf-[ ) \ u'^,n}\d\ Dog (Pf))- The 
hadronic matrix element {D*^{Pfi)\u'~fn'~f^d\D2s{Pf)) can be written as 321] 



{D*%Pf^)\uj,j,u\D,s{Pf)) 



■Tr 



Mf 



(10) 



We give the relation of wavefunctions and the strong coupling constants of different strong 
decays in appendix B. 

IV. NUMBER RESULTS AND DISCUSSIONS 



A. Semi-leptonic decays 



In order to fix the Cornell potential in equation ( lAllI) and masses of quarks, we take 
these parameters: a = e = 2.7183, A = 0.21 GeV^, Aqcd = 0.27 GeV, a = 0.06 GeV, 
rrih = 4.96 GeV, = 1.62 GeV, = 0.5 GeV, m„ = 0.305 GeV, ma = 0.311 GeV, etc 



2l| . which are best to fit the mass spectra of ground states Dq, Bq and other heavy mesons. 



In the previous letter [23|, we have obtained masses of ground states and 25' states by solving 
the BS equation. But we find that the mass of D^{2S) in 23|] is smaller than the one of 
D(2550)'^. In order to fit the experimental value of Z^(2550)°, we change the free parameter 
Vo(-D°) to Vo{D^{2S)). For other heavy-light 25* states that do not have experimental data, 
we vary Vo(lS') that has the same value as AVq = Vo{D^{2S)) — Vo{D^) to obtain the masses 
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TABLE II: Masses of the 15 and 25" states (in unit of MeV), 'Ex' means the experimental data 
from PDG 33] and 9|], 'Th' means our prediction. 



Th Ex Th Ex 

D-{1S) 1869.4 1869.6±0.16 D°{1S) 1865.0 1864.83±0.14 

D- (25*) 2560 ± 110 D^{2S) 2550 ± 109 2539.4 ± 4.5 ± 6.8 [9] 

M{2S) - M{1S) 691 ± 110 M(25) - M(15) 685 ± 109 674.6 ± 4.5 ± 6.8 

50(15') 5279.5 5279.5±0.3 D-{1S) 1968.2 1968.47±0.33 

B^{2S) 5930 ±279 Dj{2S) 2641 ± 123 

M(25) - M{1S) 651 ± 279 M(25) - M(15) 673 ± 123 

B+{1S) 5279.0 5279.17±0.29 S°(15) 5367.9 5366.3±0.6 

B+{2S) 5930 ± 280 ^°(25) 6020 ± 281 

M(25) - M{IS) 651 ± 280 M(25) - M(15) 652 ± 281 




FIG. 4: The wavefunctions of IS and the 25 of D meson with JP = 0-. 

spectra of other heavy-light 25* states. Then, we obtain the mass spectra of different mesons 
in TABLE [TTl Varying all the input parameters(A, Aqcd, ^tc) simultaneously within 5% of 
the central values, we obtain the uncertainties of masses. To show the numerical results of 
wavefunctions explicitly, we plot the wavefunctions of IS* and 25* states for D~ meson with 
= 0~ in figure |H For semi-leptonic decays, we need to input the CKM matrix elements: 
Vcb = 0.0406, Vcd = 0.23, Vcs = 0.97334 and the lifetimes of mesons: r^o = 1.53 x lO-^^^, 
ra+ = 1.638 x lO^^^^^ ^ o.453 x lO'^^^^ ^ ^ IQ-^^s, which are taken from PDG 
37|. 

In figure O as an example, we plot the form factors of the decay — )■ D~{2S)£^i'£, as 
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B%D"(1S) 






-f 





B%D"(2S) 









t -t GeV' 



t -t GeV' 



FIG. 5: The form factors of B° ^ D-{lS)i+i^i and D~{2S)i+i^e. 
TABLE III: The decay widths and branching ratios of exclusive semi-leptonic decay modes. 



Modes 



r( 10"^^ GeV) Br( 10' 



Modes 



Br( %] 



B+ D°{2S)i+Ui 2.05 ± 0.61 
B° D- {2S)i+iyi 1.95 ± 0.54 
B^ ^ D-{2S)e-^ue 4.3 ±1.2 



5.1 ±1.5 B+ D'^ilSy+ue lAr^2.2\38 
4.5 ±1.3 S° ^ L»-(15)^+t'f 1.3 ~ 2.0 [38, 
9.9 ±2.7 5° ^ L»7(lS)^+z^^ 1.4~ 1.7 [30, 



B+ B^{2S)l+vi> 0.0171 ± 0.0086 0.0120 ± 0.0060 B+ B^{lS)i+vi> 0.090 ~ 0.11 
B+ ^ B^^{2Sy+Ui 0.052 ±0.023 0.037 ± 0.016 B+ ^ B^^{lS)l+Vi 1.2 ~ 1.6 



a comparison, and we also show the form factors oi ^ D {lS)l^vi in the same method, 
where t = {P - Pff = + Mj - 2MEf and is the maximum of t. One can see that the 
values of form factors of — )■ D^{2S) are much smaller than that of B^ — 1^^(15'), and 
they have different shapes in figures. The reasons of these differences, especially the different 
shape, are mainly caused by the different wavefunctions of 15 and 25* states shown in figure 
m the numerical values of the wavefunctions for D~{1S) are positive and decrease along with 
the increased momentum while there is a node structure in D~~{2S) wavefunction, after 
the node, the wavefunctions are negative, and whose negative values increase along with 
momentum |g|. This negative part wavefunctions are responsible for the small decay width 
and special shape of form factors; similar behaviors of the form factors are also obtained in 



16 



17|. 



In TABLE IIIH we show the semi-leptonic decay widths and branching ratios with final 
mesons being 25* states, also the ones of corresponding 15* states in the same method for com- 
parison. With the same initial particle and same CKM matrix element values, the branching 
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ratios of the 25* channels are much smaller than the ones of the IS* channels. This can be 
understood by the differences of phase spaces and the node structure of wavefunctions of 2S 
state. The uncertainties of masses and decay widths shown in TABLE |TT] and TABLE IIIII 
are very large, some of them are almost 30%. The large uncertainties not only come from 
the uncertainties of phase spaces, but also from the variation of the node of the 25" wave 
function, that means a little change of node location will result in large uncertainties. 

Although compared with ground-state cases, the production ratios of 25* states are very 
small, the branching ratios of 5° D'{2S)e-^ue and B+ L'°(25)£+z/^ around lO'^ ~ 10"^ 
are very considerable. They are much larger than that of some rare decay modes and 
are accessible in the current B decay data. For the channel of — )■ Dj{2S)i'^i'£, the 
branching ratio of order 10"^ can also be accessible in the near future. But for the case 
of B^ — > B^{2S)i^i'i, due to the small Vcd and the phase space, we obtain narrow decay 
width and small branching ratio. For the decay B^ — ^ B^{2S)i^i>£, Vcs is large, while the 
phase space is very small, so the decay rate is small too. We also point out that small 
phase space and the node structure result in almost 50% uncertainties for the channels of 
B+ B%2S)i+Ui and B+ B^^{2S)£+ue. 



B. Strong decays 



These heavy-light 25* states can be detected experimentally through their strong decays. 
In the previous letter, we have calculated the strong decays of some 25* mesons |23|. But 
the predicted mass of D^{2S) is smaller than the experimental data of 1^(2550)°, which 
we have analyzed the reasons in the introduction, so we got narrow widths of Bq{2S) and 
Dq{2S) in terms of their OZI-allowed strong decays. In this paper, taking new parameters, 
we re-calculate the strong decays of Bg{2S) and Dg{2S) with large masses. For example, we 
calculate the strong decay D*^ — )■ D^n^ by the low-energy approximation, as well as the 
impulse approximation, and obtain the strong coupling constant (jr)* r)Tr= 19-8(18.3)^, which 



.3 ± 1.9 



39| and the result of 



40|. For ^ 5+7r-^hich 



is very closed to the experimental value Gd*d-k = 17.9 ± Ci. 
quenched lattice QCD calculation Gr)*Dn = 18.8 ± 2.3jl2;o 

is phase space suppressed, but as one more test of our method, similar to |33|, |4l|, |42| , we 
consider final vr as a soft-pion, and obtain the strong coupling constant 6*^.^7^=52.7(37.2)^, 



^ The first value come from tlie low-energy approximation, and the value in parentheses comes from the 
impulse approximation. 
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the impulse approximation result 37.2 is close to the same method result Gb*Bw = 30.0li;4 



m 



33|, but come with a discrepancy may come from the further instantaneous approach. 
Then, for the first radial excited states, we calculate the transition matrix elements of the 
strong decay channels, 0^(25) — t- 1^0^,0+0^, and obtain the strong coupling constants^: 

Gdo{2S)d*o^o = 3.71 ± 0.10(3.96 ± 0.34), G^o (25)^5* = 5.24 ± 0.25(5.62 ± 0.45), 
^50(25)5*0^0 = 0.82 ± 0.13(0.340 ± 0.032) GeV, ^50(25)^;-^+ = 1-10 ± 0.33(0.58 ± 0.15) GeV, 
^^-(25)^-^0 = 3.71 ± 0.10(3.91 ± 0.22), Gb~(2S)B*o^- = 5.27 ± 0.22(5.69 ± 0.60), 
G'd-(2S)d;-^o = 0.70 ± 0.12(0.40 ± 0.10) GeV, Go-(25)ds'V- = 0.92 ± 0.35(0.63 ± 0.22) GeV, 
Gd7{2S)D'0k- = 7-76 ± 0.65(6.76 ± 0.50), G'^-(25)d*-/?o = 7-77 ± 0.60(6.68 ± 0.47), 
GB+{2S)B*+nO = 9.93 ± 0.25(11.02 ± 0.31), Gb+(2s)b-o^+ = 14.04 ± 0.33(15.51 ± 0.62), 
GB+{2S)B;+nO = 1-32 ± 0.10(1.62 ± 0.12) GeV, ^^+(25)5*0^+ = 2.01 ± 0.25(2.22 ± 0.30) GeV, 
^^0(25)^.0^0 = 9.89 ± 0.30(11.02 ± 0.39) , Gbo(^2S)B'+^- = 14.07 ± 0.34(15.72 ± 0.56), 
^^0(25)^50^0 = 1.33 ± 0.11(1.62 ± 0.12) GeV, ^^0(25)5*+^- = 1-88 ± 0.22(2.34 ± 0.25) GeV, 
Gbs{2S)b*+k- = 20.09 ± 0.44(17.81 ± 0.68), ^5425)5*0^0 = 20.22 ± 0.28(17.36 ± 0.42). (11) 

Tables [IV] and |V] show the strong decay widths. Comparing equation (ITT]) with table \P7\ 
we find that the uncertainties of strong decay widths are very large, even if there are small 
uncertainties of strong coupling constants. This indicates that the predicted decay widths 
are very sensitive to the mass (or kinematic range). The two methods adopted in this 
paper, i.e. the low-energy approximation and the the impulse approximation, obtained 
similar strong decay widths, except some of the channels with P-wave involved, which have 
large discrepancies. As one can see, these P-wave-involved channels have much smaller 
phase spaces than other channels, which show that these two methods give different results 
in small phase space. And we pointed out one more time that the large uncertainties in 
table [IV] and JV] show that the decay widths are very sensitive to the node structures of 25" 
states. 

We find that there are two dominant OZI-allowed decay channels for all the heavy-light 25* 
states in table IIV] both of which are the case of 0^(25') — )■ 1^0^. Other available OZI-allowed 
strong decay channels (0~(2S') — )■ 0'''0~, 2"'"0~) include heavy P-wave meson (0^ or 2"^) in the 
final states, e.g., D^{2S) P';5(2400)°7r°, where P';;(2400)° is a scalar (0+) meson. Compared 
with the decays whose final states are all S'-wave meson (1~ and 0~), that includes P-wave 
meson in the final state and has a small decay width. In table [IV] we do not show the channels 
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TABLE IV: The decay widths (in unit of MeV) of strong decay modes. Ti and r2 come from the 
impulse approximation and the low-energy approximation, respectively. 



Modes 


Ti 


r2 


Modes 


Ti 


r2 


D°{2S) -> D*°7r° 


15.2 ±6.9 


14.1 ±7.5 


D^{2S) D*-7r+ 


31± 10 


27± 15 


5" (25') ^ I)*(2400)°7r" 


0.182 ±0.058 


0.83 ±0.30 


D°(2S) D*(2A00)-t:+ 


0.35 ±0.20 


1.39 ±0.57 


£'-(25') ^ D*-Tr° 


16.4 ±6.6 


14.5 ±7.5 


D-{2S) D*°n- 


33± 10 


30± 15 


D-{2S) i:'S(2400)-7rO 


0.21 ±0.10 


0.62 ±0.30 


D-{2S) L»S(2400)°7r- 


0.52 ±0.32 


1.10 ±0.62 


D-{2S) ^ D*°K- 


19 ±20 


25 ±26 


D-{2S) D*-K" 


17± 18 


24 ±25 


B+{2S) B*+TT° 


30.1 ±4.0 


24.1 ±6.0 


B+{2S) B*%+ 


59± 13 


48± 12 


B+{2S) B*+7r" 


0.85 ±0.25 


0.45 ±0.15 


B+{2S) B^°7r+ 


1.61 ±0.41 


1.04 ± 0.31 


B°{2S) -> B*°7r" 


30.1 ±4.8 


24.0 ±4.4 


B°{2S) B*+7r- 


60± 13 


48± 12 


B°{2S) B^^tt" 


0.85 ±0.24 


0.46 ±0.15 


B"{2S) B;+tt- 


1.81 ±0.50 


0.91 ±0.31 


Bsi2S) B*+K- 


42± 18 


55± 19 


B,{2S) B*°K" 


41± 18 


54± 19 



that involve heavy tensor meson {J^ = 2+), which has an ignorable decay width caused by 
the extremely narrow kinematic range. There are not enough experimental data for the 
scalars mesons (J^ = 0^). For example, Z)o(2400)° has been confirmed, but Dq(2400)~ is 
missing. So we choose the same mass for Z)g° and Dq~ , i.e. Mf)*(^24oo)'> = 2.296 ±0.095 GeV, 
and the same mass for Bq^ and Bq^ , i.e. M^*o = 5.660 ± 0.266 GeV. For the decays of 
Ds{2S) and Bs{2S), the channels involving P-wave state in the final states are kinematic 
forbidden. 

Since we have obtained higher masses for heavy-light 2S states than the ones of the 
previous letter 23j, we obtain much broader widths, but the estimated full width of D^{2S) 
is still narrower than the experimental data of D{2550)^. Our prediction is Ti)0(^2S) = 
47 ± 17(43 ± 23) MeV, while from data is 130 ± 12 ± 13 MeV [9;]. The new state ^(2550)° 
should be the state D^{2S) around the 2540 MeV, since other states can be ruled out by 
the masses or decay modes. For example, P-wave D*j states whose masses are closed to 
2540 MeV have strong decay productions including K or Ds but not D*7i. The P-wave 
Dj states, which have masses about 100 MeV lower than 2540 MeV, have strong decays 
different from 0~(2S'). For Z^-mesons, which are 1^{2S), 1^(2D), or the mixing of 1~(2S') 
and 1^{2D), their masses should be higher than 2540 MeV and we already have the candidate 
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TABLE V: The estimated full widths (in unit of MeV) of the first radial excited states. The first 
and second value (in brackets) come from low-energy approximation and impulse approximation, 
respectively. 



D°{2S) D-{2S) D-{2S) B°{2S) B+{2S) B,{2S) 

43 ± 23(47 ± 17) 46 ± 23(50 ± 17) 49 ± 51(36 ± 38) 74 ± 17(91 ± 18) 73 ± 19(93 ± 19) 109 ± 38(83 ± 36) 

D*(2600) [9] that have more strong decay channels: 0~1~, 0~0~, 1~0^. Therefore, 

the new state D(2550)'' is the state D'^{2S). The discrepancy between the theoretical and 
experimental result comes from that we only considered the dominant OZI-allowed strong 
decays. Finally, we should point out that the experimental data may vary along with new 
more precise detections. 

C. Product of semi-leptonic decay ratio and cascaded strong decay ratio 

We have calculated Bg/ semi-leptonic decay to Bq(2S) and Dg(2S), where all the states 
are on mass shells; we also calculated the main strong decays of Bq{2S) and Dq{2S), e.g., 
we obtain the branching ratio of the semi-leptonic decay Br{B^ — t- D^{2S)l^Vi) = (5.1 ± 
1.5) X 10-^ the strong decay branching ratio Br{D^{2S) D*-n+) = ^^"'^^'j^^^p"" = 

0.63±0.48(0.66±0.15) and Br{D^{2S) D*%°) = ^"^"r^p^J*""" = 0.33±0.24(0.32±0.11). 
To reduce the influence of the discrepancies caused by the theoretical strong decay widths, 
and the ground IS* states D^*\ B^*\ and their decays that are well known in the experiment, 
we multiply the branching ratio of the semi-leptonic decay with the strong decay branching 
ratio, which show the ability of experiment to detect the missing 25" states, but we ignore 
the reconstructed efficiencies of events in experiment, and the products of ratios are 

Br{B^ D^{2S)tvf) x 5r(D°(2S) ^ D*"7r+) ^ [3.2 ± 2.6(3.4 ± 1.9)] x 10"^ 
Br{B+ D%2S)i^ug) x Br{D\2S) D*V) ^ [1.7 ± 1.4(1.6 ± 1.0)] x 10'\ (12) 

Br{B^ D-{2S)tvi) x Br{D-{2S) D*\-) ^ [3.0 ± 2.3(3.0 ± 1.6)] x 10-^ 
Br{B^ D-{2S)tvi) x Br{D-{2S) D*-n^) ^ [1.4 ± 1.1(1.5 ± 0.9)] x 10"^ (13) 

Br{B° D-{2S)tvi) x Br{Dj{2S) D*°K-) ^ [5.1 ± 7.6(5.2 ± 7.9)] x 10-^ 
Er(5° ^ D;{2S)e+ue) X Br{D;{2S) D*- K°) ^ [4.8 ± 7.3(4.7 ± 7.1)] x 10~^ (14) 
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Sr(S+ ^ B°{2S)tue) x Br{B^{2S) B*+n-) ^ [0.79 ± 0.47(0.79 ± 0.46)] x 10'^, 
Br{B+ B^{2S)tvi) x Br{B\2S) 5*V) ^ [0.39 ± 0.23(0.40 ± 0.22)] x 10-^(15) 

Br{B+ Bl{2S)tvi) x Br{Bs{2S) -> B*+K-) ^ [1.9 ± 1.5(1.9 ± 1.4)] x 10"^ 
Br{B+ B^^{2S)i+ue) x 5r(5,(25) ^ 5*°/?°) ^ [1.8 ± 1.5(1.8 ± 1.4)] x 10"^ (16) 

The decays of B^ and B^ have ratios of order lO"'*, which can be analyzed with current 
data at B-factories. The decay of Bs that has ratios of order 10"'' may be observed in the 
future, while the Be decay ratio of order 10~^ is hard to reach experimentally. 

In summary, we have studied the productions of Dq{2S) and Bq{2S) in the exclusive 
semi-leptonic Bg' decays and the strong decays of Dq{2S) and Bq{2S). Some of these decays 
have the branching ratios of order 10~^, which could be measured currently in experiments. 
For examples, the ratios Br{B+ D\2S)i+ue) x Br{D\2S) D*tx) ^ [4.9 ± 4.0(5.0 ± 
2.9)] X 10"^ and Sr(S° D-{2S)i+ue) x Br{D-{2S) D*n) ^ [4.4±3.4(4.5±2.5)] x 10"^ 
are relatively large in B decays, which could be detected by the two current 5-factories. For 
Ds{2S), the ratio Br{B^^ D;{2S)£+iye)x Br{D;{2S) D*K) ^ [9.9± 14.9(9.9± 15.0)] x 
10~^ is also not small, which will be reached in the future. We have also given the strong 
coupling constants of Dq(2S) and Bq{2S), which maybe observed experimentally. Although 
similar to other models' results, our calculation also gave a smaller full decay width than 
the experimental data, but the recent detected -D(2550)'^ is very likely the D^{2S) state. 
Acknowledgements This work was supported in part by the National Natural 
Science Foundation of China (NSFC) under Grant numbers 10875032 and 11175051. 



Appendix A: Instantaneous BS Equation 

In this section, we briefly review the BS equation and its instantaneous one, i.e., the 
Salpeter equation. 



The BS equation is read as [43 1 



r d^k 

(^1 - mi)x(g)(^2 + m2) = ij J^ViP^ q)x{k) , (Al) 

where x(q') is the BS wavefunction, V{P, k, q) is the interaction kernel between the quark 
and anti-quark, and piandp2 are the momenta of the quark 1 and anti-quark 2. 
We divide the relative momentum q into two parts, i.e., gy and q±: 
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In the instantaneous approach, the kernel V{P, k, q) takes the simple form 4J]: 



V{P,k,q)^V{\k-q\). 

Let us introduce the notations fp{qj_) and 77(9^^) for a three-dimensional wavefunction as 
follows: 

dqp^ 



ri{q[ 



dk I 



(2vr)^ 



:V'(A:±,g±)</?p(A;^^) . 



Then, the BS equation can be rewritten as 



The propagators of the two constituents can be decomposed as: 



SiiPi 



+ 



with 



J{i)qp + oiiM — uji + ie J{i)qp + a^M + Ui — ie 



— ±J{i){mi + 



where z = 1, 2 for quark and anti-quark, respectively, and J{i) = ( 
We introduce the notations as 

ipi^{q±) = A%{q±)j^Ml±)j^^%i'l±) ■ 
With contour integration over qp on both sides of equation flASp . we obtain 

. Aipi(l±)Vp{(l±)Mpi(l±) Arp(g±)r7p(g±)A2p(g±) 



(M — cui — Ci;2) 



(M + + 002) 



and the full Salpeter equation 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



{M-Ui- uj2)^+^{q±) = Atp{q±)Vp{q±)K{q±) , 
{M + Ui+U2)y^;'{q±) = -Aip{q±)Vp{q±)A2p{q±) , 

¥>t'(l^) = Vp^(^^) = ^- (A7) 
For the different J^*" (or J^) states, we give the general form of wavefunctions (we will 
talk about them in appendix B). Reduce the wavefunctions by the last equation of (1A7p . 
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and then solve the first and second equations of (1A7P to obtain the wavefunctions and mass 
spectrum. We have discussed the solution of the Salpeter equation in detail in {21, 29 1. 
The normalization condition for the BS wavefunction is 



27r2 



^ M ^ M 



2Pn . 



(A8) 



In our model, the instantaneous interaction kernel V is the Cornell potential, which is 
the sum of a linear scalar interaction and a vector interaction: 



4 a 

y{r) = K(r) + Vo + 7o ® 7°K(r) = Xr + Vo - 7o ® 7^ — 

6 r 



(A9) 



where A is the string constant and as{q) is the running coupling constant. In order to fit 
the data of heavy quarkonia, a constant Vq is often added to the confine potential. One can 
see that V^{r) diverges at r = 0; we introduce a factor e~"^' to avoid the divergence 



VJr) 



A 



a 



(1 - e-"0 , K(r) = - 



4 a, 
3 r 



(AlO) 



It is easy to know that when ar <^ 1, the potential becomes equation flA9l) . In the 
momentum space and the center of mass system of the bound state, the potential reads : 



A 



Vsiq) = -(- + Vo)6%q) 



A 1 . 2 



r2 r;^ 



where the running coupling constant as{q) is 

127r 



as{q) 



37r2 + a'') 



(All) 



33 - 2Nf log(a + ^ 



^QCD 



We introduce a small parameter a to avoid the divergence in the denominator. The constants 
A, a, Vo and Aqcd are the parameters that characterize the potential. Nf = 3 for bq (and 
cq) system. 



Appendix B: Wavefunctions for different states 

We know that form factors of semi-leptonic decay and strong coupling constants of strong 
decay are related to wavefunctions of different states in section II and III. In this section, 
we give the wavefunctions of the different states and obtain the form factors and strong 
coupling constants. 
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a). For Bqi meson with quantum numbers = 

The general form for the relativistic wavefunction of pseudoscalar meson can be written 
as [291 

^^1 



M 



75, 



(Bl) 



where M is the mass of the pseudoscalar meson, and fi{q) are functions of \q\^. Due to the 
last two equations of (1A7P : (p^Z = ^Pq^ = 0, we have 



m2UJi + miuj2 



fM 



m2UJi + miW2 



(B2) 



where mi and 1712 and Ui = ymf + (p, 002 = ym^ + are the masses and the energies of 
quark and anti-quark in Bg' mesons, = Q — ' P/M^)P, and = — |g|^. Then, there 
are only two independent unknown wavefunctions /i(g) and /2(g) in equation fIBip : 

M{ui — UJ2) 



1712001 + miUJ2 



OJi + U2 



75- 



(B3) 



m2UJi + miuj2- 

The numerical values of radial wavefunctions /i and /2 and eigenvalue M can be obtained 
by solving the first two Salpeter equations in equation ( IA7p . 

According to equation flA6p . the relativistic positive wavefunction of pseudoscalar 0" 
state in center of mass system can be written as 



29| 



O2 H h O'^ I + Oa 



75, 



(B4) 



where the 6jS (i = 1, 2, 3, 4) are related to the original radial wavefunction /j, quark masses 
rrii, quark energy Wi {i = 1, 2) and meson mass M: 



bl = ^ [ fM + f2iq} 



nil + 1^2 
UJi + 002 



Wi + 002 
nil + ^2 



(mi - m2) 
miW2 + m2ix'i 



[UJl + U2) 



{miuj2 + m2Cc;i) ' 



b). For Bq{2S) and Dq{2S) mesons with quantum numbers = 0~ 

Because the 2S state mesons have the same quantum numbers as Bq/, the wavefunction 
of 25* state mesons are similar to equation flB4p . 



ai 



02 



It 



as 



04- 



Mf 



75, 



Mf 



ai 



m'l + m'g 



02 



/ ' 



m'l + m2 



as 



{m'l — ni'2) 

— ; — ~, — r~r'"4 

miCi;2 + m2Ci;i 



(B5) 
{uj'i+u'2) 



[ni'iu'2 + ni'2Uj'i, 
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where Mf, Pf and fHqj) are the mass, momentum and the radial wavefunction of 25* state 
mesons, respectively, m'^^m^ and uj'i - 



im'i + qj,uj2 = + are the masses and the 

energies of quark and anti-quark in 25* state mesons, respectively. 

According to the equations ([9]), flB4p and (IBSp . the form factors of Bq' semi-leptonic 
decays to 25' state mesons can be written as 



1 f (fq iaibi 
2J (27r)3 MM 



/ 



(02^2^/ + b^q ■ r cos 6 — 0462-^^0'^ — a^b^i^Efq ■ r cos 6)a[ 



+a2,Mf{h2,q^ + -E/a'i — h^q ■ rcos6a\) + M(l + a^b^q^ + aib2Efa[ — a2,Mfa[] 
-M(l - ^)(a4&2^/ - hEf - asMf + a2hiMf + 0464? ■ rcos6') 



(B6) 



/- = - y J27iy MAdf [('^2&2M/ + fegg ■ r COS 61 - aib2Eja[ - a^b^Efq -r cos 9)a[ 
+a3Mf{b^q^ + Efa'i — b^q ■ rcos6a\) — M(l + a^b^q'^ + aib2Efa[ — asM^a']^) 
+M(1 + ^)(a4&2-E/ - - asM/ + a2biMf + 0464? ■ r cos6 



(B7) 



where -E/ = ^J Mj + , g ■ r = |(f ■ r|, ^ is the angle between q and r. 
c). For B* and mesons with quantum numbers = 1 

The relativistic positive wavefunction of 1^ state can be written as 



■ e^^^(c5 + C6 /'/I + C7 + C8 f/l), (B8) 

where we first defined the parameter rii that is functions of /" (1^ wave functions): 
ni = h [Qfi) - /e ilnJjZJ^, ^> ^2 = /s (g/i) - /e W/ 



K' + m'2'; 
ns = /^'(g/i) + mqn 



+^2) 



Then, we defined the parameters q that are functions of /" and n,: 



M/i 1 (m'/ + m'2') 

ci = — ^^i,C2 = --^-n-. — —ni^c^ = ^ , ; , — ^r^ni,C4 = 



C5 



C7 



2 (a;'/ + uj'i) 



2 (m'iu!^ + m'2'cj'/) ' 
1 



+ Tn'2 / 71 j-2 2 \ ^1 — ^2 / 71 /r2 2 \ 



2M/1 (wi^a + m'{m2 + qj 

nmMfi 

2M/1 (m'M + m'2'a;l') ' ^ 2M|i m'/ + m'2' 



1 u'l + w^' 



77^3-/5(9/1 



w'/ + w^' 



(m'l' + m'2')(a;'M + m'(m'i - q),^] 
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According to equations (fTOl) . (IBSP and (IBSp . the strong coupling constant for (25") — )■ 
1^0^ can be written as 



1 /" d'^Q f 1 

G'o-{25Hi-o- ^ J (2^^^M ilF ^^^^^^f + a3C4Mji)|g>p - aiC2Mfqf ■ Pji 
+a.MQf ■ Pfi? - Ci{Mf + a^qf ■ Pfi) + c^E f^al^E f^M f - Mj, ) + a^c^\qf\^Ef^al{Mf - Ej 



aiCs\qf\^Efia2{EfiMf - Mj^) + {cjEfi - a^c^Eji + a^c^Ej^ + a^csMj + a3C6EfiMf)qf ■ Pfia^ 
-a^CsEfia'^iqf ■ Pj^f + Cr^/iaf |P^ip + agCrg/ ■ PfiafEf^{Mf - Ef{) 



+a4C8S;ia2'"(M/ - ^/i) + a2{c2Mj^ + C3(g/ ■ P/i + iP/ira'a. 

+Ef,al{c,{Ef, - Mj) + csEj,{qf ■ Pf, + Ej.a'i - M],a'^))) ] - ^^iM [a,c5(M^ - Ef,) 

Mf\Pfi\ 

+ (a4C2 + C6)(M^i - EfiMf) + a4C8|g/|^(M^i - EfiMf) + (0403 + 0405 - C7 - a2C^Efi)qf ■ Pji 
+a4Cs\qf ■ P/iP + (0403 - C7)Ej^a2 - aaCg^Jia'a + (C7 - aiC3)Mj^a2 + aaCg^iMjia's 
+a4C8g> ■ Pfia'iEfi{Ef^ - Mf^) + a3(ci(My - Eyi) + (04^/1 - CgMy).?/ ■ P^i 
+c,Ef,al\Pf,\'' + C7(i?/i - My)(|g->|2 + g) ■ P^ia^'))] } (B9) 



where fx {X is vr or i^') is the decay constant for the hght meson, a'j — ^„^^„ . 
d). For P* and P* mesons with quantum numbers = 0^ 
The relativistic positive energy wavefunction of 0''" can be written as 

^oV(g7l) = C,Un^ + C2^jJ^ + C3 + C4^), (BIO) 

where the parameters Cj are the functions of /" and /2 (0^ wavefunction), which are defined 

as 

Using equations ffTU]) . fIBSp and fIBlOp . the strong couphng constant for 0^(25) — )■ 0+0~ 
can be written as 

Go-{25Ho+o- =J^I (2 J)^3 ~M^'r {(a4C2P/iM/ - a^Ef^Mji + a^MfMfi - a^C2M]^)\qff 
+ {Mfi - a^C3Mf,)qf ■ P^i - a4C2|g/ ■ P/i|' + c^iEf.Mf - Mj, + a^Mfqj ■ Pf{) 
+ \Pfi\^Mfia2 + (a4C2P/i + Mfia3){Mf - Efi)qf ■ Pjia'^ + a2{c3Mfi{Efi - Mf) 
+C2P/i(g>-P;i + |P;i|V2)} (Bll) 
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